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ABSTRACT 

The power spectrum (PS) of mass density fluctuations, independent of "bi- 
asing", is estimated from the Mark III Catalog of Peculiar Velocities of galaxies 
using Bayesian statistics. A parametric model is assumed for the PS, and the free 
parameters are determined by maximizing the posterior probability of the model 
given the data. The method has been tested using detailed mock catalogs. We 
use generalized CDM models with and without COBE normalization. 

The robust result for all the models tested is P(k)Q} 2 = (4.1 ± 0.7) X 
10 3 (h~ 1 M.pc) 3 at k = 0.1/iMpc -1 , with the peak constrained to the range 
0.03 < k < 0.06 h Mpc -1 . It is consistent with a direct computation of the PS 
(Kolatt & Dekel 1996). When compared to galaxy- density surveys, the implied 
values for (3 (= ft°' 6 /&) are of order unity to within 25%. 

A T-shape model, free of COBE normalization, is constrained by the velocity 
data toT = 0.5 ± 0.15 and a 8 tt - 6 = 0.85 ± 0.1. Within the family of COBE- 
normalized CDM models, the best tilted model (0 = 1, 0.84^ 65 ) and the 
best open model (n = 1, ft 0.75h^ ' ) are more likely than the best A model 
(n = 1, A = 1 - ft, ft « 0.70/?," 1 - 2 ). The most likely CDM model with ft < 1 is 
found to be of ft = 1, h 0.75, and a tilted spectrum of n = 0.8 ± 0.02 with 
tensor fluctuations. The tightest constraint for the tilted-A models is of the sort 
Qhl^n^ = 0.7 ± 0.08, with v = 3.8 and 1.85 with and without tensor fluctuations 
respectively. 

Subject headings: cosmology: theory — cosmology: observation — dark matter 
— galaxies: clustering — galaxies: distances and redshifts — large scale structure 
of universe 



1. INTRODUCTION 

In the standard picture of cosmology, the structure on large scales originated from 
small-amplitude initial density fluctuations that were amplified by gravitational instability. 
These initial fluctuations are assumed to be a Gaussian random field, solely characterized 
by it's power spectrum. On large scales, the fluctuations are linear even at late times, so 
that the power spectrum preserves it's original shape. This makes it a very useful statistics 
for large-scale structure. 

The power spectra of galaxy density were derived for many different samples, in two 
angular dimensions or in three dimensions from redshift space. Unfortunately, these power 
spectra correspond to objects that are not necessarily unbiased tracers of the underlying 
mass distribution, and it is the mass distribution that is directly related to theory (e.g. 
Dekel & Rees 1987 for a review on "galaxy biasing"). Clear evidence for this bias is 
provided by the fact that galaxies of different types are observed to cluster differently 
(e.g. Dressier 1980). It would therefore be naive to assume that any of the galaxy power- 
spectra directly reflects the mass PS. Furthermore, in estimates of the galaxy PS from 
redshift surveys, uncertainties also arise when correcting for redshift distortions (Kaiser 
1987, Zaroubi and Hoffman 1996). For these reasons, one wishes to measure the mass PS 
directly from dynamical data, bypassing the complex galaxy-biasing issues and the need to 
correct for redshift distortions. In principle, such dynamical information can be provided 
by peculiar velocities, by gravitational lensing effects, or by fluctuations in the cosmic 
microwave background (CMB). In particular, the accumulating catalogs of galaxy peculiar 
velocities enable a direct determination of the mass PS under the natural assumption that 
the galaxies are unbiased tracers of the large-scale, gravitationally-induced velocity field. 

Here, the PS is computed from the Mark III Catalog of Peculiar Velocities (Willick et 
al. 1995 WI; 1996a WII; 1996b Will), which consists of more than 3000 galaxies. It was 
compiled from several different data sets of spiral and elliptical/SO galaxies with distances 
inferred by the forward Tully- Fisher and D n — a methods, which were re-calibrated and 
self-consistently put together as a homogeneous catalog for velocity analysis. The catalog 
provides radial peculiar velocities and inferred distances, all properly corrected for inho- 
mogeneous Malmquist bias, for ~ 1200 objects, ranging from isolated field galaxies to rich 
clusters. The associated errors are on the order of 17 — 21% of the distance per galaxy. 
The sparse and inhomogeneous sampling is another source of error. 

These data allow a reasonable recovery of the dynamical fields with ~ 12/?. _1 Mpc 
smoothing in a sphere of radius ~ 60/?. _1 Mpc about the Local Group, extending to ~ 
80/i -1 Mpc in certain regions. The POTENT method (Bertschinger & Dekel 1989; Dekel, 
Bertschinger & Faber 1990; Dekel 1994) attempts a recovery of the underlying density field 
with fixed Gaussian smoothing within this volume. In an associated paper, Kolatt and 
Dekel (1996, KD) have computed the mass PS from the smoothed density field recovered 
by POTENT from Mark III. The limitations of the data introduce severe systematic errors, 
that were modeled via Monte-Carlo mock catalogs and then used to correct the measured 
PS. Since the KD results naturally involve uncertainties, an independent estimate of the 
PS, using a very different method, is useful. 

Our purpose is to estimate the mass PS directly from the peculiar velocities of the 



2 



Mark III catalog, by means of a likelihood analysis. The non-local nature of the peculiar 
velocities, i.e. being influenced by the mass distribution in a whole neighborhood, allows 
one to probe scales somewhat larger than those probed by the density field. For example, 
the effect of a bulk velocity across the entire volume is not evident if only the density field is 
considered. For a similar reason, the velocity field is expected to obey linear theory better 
than the density field smoothed on a comparable scale, and so to be closer to a Gaussian 
field. Our approach here does not involve any explicit window function, weighting or 
smoothing, nor does it require artificial binning of the PS. In addition, it automatically 
underweights noisy, unreliable data. 

The data analyzed here are especially suited for Bayesian analysis. The sparse and 
inhomogeneous sampling of a random Gaussian field with Gaussian errors yields a multi- 
variate Gaussian data set. The corresponding posterior probability distribution function 
(PDF) is a multivariate Gaussian that is completely determined by the assumed PS and 
the assumed covariance matrix of errors. Under these conditions one can write the joint 
PDF of the model PS and the underlying velocity or density field, and then simultaneously 
estimate the PS model parameters and recover the "Wiener filter" solution of the fields 
(Zaroubi et al. 1995). In an associated paper (Zaroubi, Hoffman & Dekel 1996), we present 
the high-resolution fields recovered from this same data set using the PS derived here. 

To apply our method, the simplifying assumptions that have to be made are that the 
peculiar velocities are drawn from a Gaussian field, that their correlation function can be 
derived from the density PS using linear theory, and that the errors are Gaussian. The 
need to assume a parametric functional form for the PS is also a limitation; one can try to 
achieve flexibility by using a large number of parameters and a variety of functional forms, 
but at the risk of sometimes making the likelihood analysis unstable (§5). 

The method is described in §2, in which the relation between the velocity correlation 
functions and the PS, and the general likelihood-analysis algorithm for computing the PS, 
are specified. The method is tested using a mock catalog in §3. The resultant power spectra 
are presented in §4, as derived from the Mark III data alone, and for generalized CDM 
models imposing COBE normalization. The associated constraints on the cosmological 
parameters are analyzed. Our conclusions are summarized and discussed in §5. 



The computation of the matter power spectrum from the peculiar velocity data by 
means of likelihood analysis requires a relation between the velocity correlation function 
and the power spectrum. Define the two-point velocity correlation (3 X 3) tensor by the 
average over all pairs of points r 4 - and rj that are separated by r = rj — 



where v^(ri) is the //. component of the peculiar velocity at r 8 -. In linear theory, it can be ex- 
pressed in terms of two scalar functions of r = |r| (Gorski 1988), parallel and perpendicular 
to the separation r, 



2. METHOD 



2.1. Velocity Correlations 



(1) 



^(r) = V±(r)8 llv + [tf ||(r) - *±(r)]r^r 



(2) 
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The spectral representation of these radial correlation functions is 



^'" (r) = ^0 }1 l^P(k)K ± , ll (kr)dk, (3) 

where K±(x) = j\{x) / x andif||(a;) = jo — 2ji(x)/x, with ji(x) the spherical Bessel function 
of order I. The cosmological dependence enters as usual in linear theory via f(£l) ~ fi ' 6 , 
and Ho is the Hubble constant. A parametric functional form of P{k) thus translates to a 
parametric form of V^j,. 

2.2. Likelihood Analysis 

Let m be the vector of model parameters and d the vector of N data points. Then 
Bayes' theorem states that the posterior probability density of a model given the data is 

, P(m)P(d|m) 

p(m|d)= p(d) (4) 

The denominator is merely a normalization constant. The probability density of the model 
parameters, P(m), is unknown, and in the absence of any other information we assume it 
is uniform within a certain range. The conditional probability of the data given the model, 
P(d|m), is the likelihood function, £(d|m). The objective in this approach, which is to 
find the set of parameters that maximizes the probability of the model given the data, is 
thus equivalent to maximizing the likelihood of the data given the model (c/. Zaroubi et 
al. 1994; Jaffe & Kaiser 1994) 

The Bayesian analysis measures the relative likelihood of different models. An absolute 
frequentist's measure of goodness of fit could be provided by the Chi-squared per degree 
of freedom, which we use as a check on the best parameters obtained by the likelihood 
analysis. 

Assuming that the velocities are a Gaussian random field, the two-point velocity 
correlation tensor\t fully characterizes the statistics of the velocity field. Define the radial- 
velocity correlation (N X N) matrix Uij by Uij = rj\trj, where i and j refer to the data 
points. Let the inferred radial peculiar velocity at r 8 - be Ui, with the corresponding error e 8 - 
also assumed to be a Gaussian random variable. The observed correlation matrix is then 
Uij = U{j + e?6ij, and the likelihood of the N data points is 



C = [(27r) JV det([/ u -)]- 1/2 exp --TutU^uj . (5) 




Given that the correlation matrix, Uij, is symmetric and positive definite, we can 
use the Cholesky decomposition method (Press et al. 1992) for computing the likelihood 
function (Eq. 5). The significant contribution of the errors to the diagonal terms makes the 
matrix especially well- conditioned for decomposition. The calculation for a given choice of 
parameters and N ~ 1200 data points takes several minutes on a Dec-Alpha workstation 
(of SpecFP92 ~ 150). 
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The likelihood function of Eq. 5 is the posterior PDF of the parameters m. It is 
also a x 2 distribution (with N degrees of freedom) with respect to the N data points, but 
it is not necessarily a % 2 distribution with respect to the parameters, and it is therefore 
difficult to assign accurate confidence levels to the parameters. This requires elaborate 
integrations over the volume encompassed by the equal- likelihood surfaces in parameter 
space. In the present paper we limit ourselves to a rough estimate of confidence levels by 
crudely approximating the PDF as a % 2 distribution in parameter space. 

Note that the quantity that can be derived from peculiar- velocity data via the linear 
approximation is f 2 (Q,)P(k), where P is the mass density PS (see Eq. 3). 

3. TESTING THE METHOD 

The Mark III catalog (WI; WII; Will) provides inferred forward TF velocities to 
~ 1200 objects, grouped from more than 3000 galaxies in order to reduce the distance un- 
certainty per object and thus reduce Malmquist bias. The velocities were further corrected 
for homogeneous and inhomogeneous Malmquist bias {e.g. Dekel 1994). 

The grouping serves us here also as a means of smoothing over nonlinear velocities. 
The PS on large-scales should not be too affected by this grouping because the individual 
galaxies in clusters enter the likelihood with low weights anyway. This is because the 
randomness of the velocities within clusters makes the noise term e 2 8ij dominate over the 
signal Uij in the observed correlation matrix Uij in Eq. 5. 

Careful testing of the method with realistic mock catalogs is essential in view of the 
large distance errors, the sparse and non-uniform sampling, the bias-correction procedures, 
and the possible non-linear and non-Gaussian effects. 

The mock Mark III catalogs are described in Kolatt et al. (1996). They are based on 
simulations whose initial conditions were extracted from a reconstruction of the smoothed 
real-space density field from the IRAS 1.2Jy redshift survey, taken back into the linear 
regime. Small-scale perturbations were added by means of constrained random realizations. 
The system was then evolved forward in time using an N-body simulation assuming = 1, 
and stopped at two alternative times, when the rms density fluctuation in a top-hat sphere 
of radius 8/?. _1 Mpc reached ag = 0.7 and later when sigmag = 1.12. The "galaxies" 
in the simulation were identified via a linear biasing scheme (b=1.35), and they were 
divided into 'spirals' and 'ellipticals' according to Dressler's morphology-density relation. 
The galaxies were assigned TF quantities (internal velocities and absolute magnitudes) 
that were Gaussianly scattered about an assumed TF relation, and were then "observed" 
following the selection criteria of the actual data sets that compose the Mark III catalog. 
The mock catalogs were grouped and corrected for biases just like the real catalog. 

The true PS of the mass in the simulation is well approximated by the functional form 

with A = 4.68 and 12.28 x 10 4 ( fe -1 Mpc) 4 , B = 8.3 and 9.25fe -1 Mpc, and a = 3.2 and 
2.8, for the ag = 0.7 and 1.12 cases respectively. 
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Figure la: Contour map of log likelihood in the a— A plane for a mock catalog based on 

the parametric model of Eq. 6 with B = 9.25 /?. _1 Mpc. Contour spacing is A[ln(£)] = 1. 

The best-fit point is marked. 



Figure lb: The true power spectrum of the simulation (heavy solid), compared with 
the best-fit solution (heavy dashed), and two power spectra whose parameters lie on the 
innermost closed contour of Fig. la. The values of ag for this contour are 1 — 1.2. 

The likelihood analysis has been applied to the mock catalogs using the parametric 
functional form of Eq. (6) as a prior. To save effort we always kept one of the parameters 
fixed and allowed only the other two to vary. Figure la shows for example a contour 
map of ln(£) for one of the ag = 1.12 mock catalogs, spanning the a- A plane with B = 
9.25 /i -1 Mpc. The contours are separated by Aln£ = 1. Maximum likelihood is obtained 
at A = Ao and a = 2.75 (compared to 2.80). Assuming a Chi-square distribution with 
two degrees of freedom, the la contour of the likelihood around the best-fit parameters is 
at ln(£) 2.1. We conclude that the recovery of the PS is excellent. 

Figure lb shows the recovered PS in comparison with the true PS of the simulation. 
They almost coincide over the whole range of scales, showing slight deviations only on very 
small scales. To illustrate the level of uncertainty, we plot for comparison several other 
power spectra that were obtained with parameter pairs that lie on the innermost contour 
about the maximum in Fig. la. It shows that the amplitude near the peak can be off by 
about 25%, and that the recovery becomes more robust at moderately smaller scales. The 
success of the recovery is similar when the other pairs of parameters are allowed to vary. 

4. RESULTS 

4.1. TheY Model 

To recover the PS from the velocity data independent of the COBE normalization, we 
use as a parametric prior the so-called T model (e.g. Efstathiou, Bond and White 1992), 

P(k) = AkT 2 (k), T(k) = (l + [ak/T + (bk/Tf/ 2 +(ck/Tf] v y 1/ \ (7) 
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with a = 6.4 h~ 1 M.pc, b = 3.0 h~ 1 M.pc, c = 1.7 h~ 1 M.pc and = 1.13. The free parameters 
to be determined by the likelihood analysis are the normalization factor A and the T 
parameter. In the context of the CDM cosmological model, T has a specific cosmological 
interpretation, T = Qh. Here, however, independently of CDM, Eq. 7 serves as a generic 
function with logarithmic slopes n = 1 and —3 on large and small scales respectively, and 
with a turnover at some intermediate wavenumber that is determined by the single shape 
parameter T. 

Figure 2a shows the contour map of In C in the T-A plane. The maximum likelihood 
values are V = 0.5 and A = 9.4 X 10 4 ( h~ 1 M.pc) i . The corresponding value of ergft ' 6 is 
0.85 ± 0.1. Figure 2b shows the best-fit PS (solid). To illustrate the uncertainty in the 
PS we also show the power spectra of five other parameter pairs that lie on the innermost 
likelihood contour about the best fit. 
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Figure 2a: Contour map of log likelihood for the T model. Contour spacing is A[ln(£)] = 
1. A in units of A = 6.28 x 10 5 ( h^Mpc) 4 . 

Figure 2b: The most likely T-model power spectrum (solid), and five other models whose 
parameters lie on the innermost contour of Fig 2a. 



4.2. COBE-Normalized CDM Models 

We now restrict our attention to the generalized family of CDM cosmological models, 
allowing variations in the cosmological parameters fi, A and h, as well as the large-scale 
PS slope n and the contribution of tensor fluctuations. Furthermore, we now impose 
the normalization implied by the two-year COBE DMR data as an additional external 
constraint. The general form of the PS in these models is 

P(k) = A C oBE(n, ft, A) T 2 (ft, n B ,h- k) k n , (8) 

where we adopt the CDM transfer function proposed by Sugiyama (1995), 

T(k) = ln(1 2 + 4 2 - 3g)) [1 + 3.89g + (W.lqf + (5.46g) 3 + (6.71 g ) 4 ] " 1/4 , (9a) 
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(9b) 

The parameters are varied, two at a time, such that they span the range of currently 
popular CDM models, including Tilted-A CDM (ft + A = 1, ft < 1, n < 1) and Tilted- 
Open CDM (A = 0, ft < 1, n < 1). We allow the possibility of nonzero tensor fluctuations 
of T/S = 7(1 — n), where the ratio is of the C 2 quadrupole moments of temperature 
angular fluctuations of the tensor and scalar modes (e.g. Turner 1993; Crittenden et al. 
1993). In all cases, the baryonic density is assumed to be ft& = 0.0125/?. -2 , which is the 
value currently favored by primordial nucleosynthesis analysis (e.g. Walker et al. 1991). 

The COBE normalization for each model has been calculated by various authors 
(Gorski et al. 1995; Sugiyama 1995; White & Bunn 1995), using different Boltzmann codes, 
different statistical analyses, and sometimes even different temperature maps. We have 
arbitrarily adopted Sugiyama's normalization as a backbone, and for models not studied 
by him we use the other results after matching them to Sugiyama's using the models that 
they have investigated in common. 

In particular, the COBE normalization is modeled by Acobe = (ft) A2 (n). For 
Tilted-A CDM models we use the fits: 

Ai(ft) = dex(7.93 - 8.33ft + 21.31ft 2 - 29.67ft 3 + 10.65ft 4 + 

15.42ft 5 - 6.04ft 6 + 13.97ft 7 + 8.61ft 8 ), (10a) 

Mn) = \f K l- 2 ;T 2 :lT\ <i<») 

A ; \ dex(-4.54 + 4.54n) (T / 0). v ; 

Thes fits are for h = 0.5, but the h dependence in the range of interest is weak, and we 
ignore it here. 

For Tilted-Open CDM model with T = the fit is: 

Ai(ft) = dex(5.75 - 1.68ft + 4.53ft 2 - 7.57ft 3 + 7.53ft 4 + 3.15ft 5 - 0.23ft 6 ), (11a) 

A 2 (n) = dex(-2.71 + 2.71™). (116) 

4.2.1 Low Omega Models 

Figure 3a shows the likelihood contour map in the ft — h plane, for the ACDM family 
of models with n = 1 (normalization by Sugiyama). The most probable parameters in 
this case (in the range ft < 1) are ft = 1 and h = 0.4. Figure 3b shows the similar 
likelihood map for OCDM with n = 1. The most probable values here are ft = 0.46 and 
h = 0.9. It is clear from the elongated contour maps that ft and h are not constrained 
very effectively independently of each other. It is a degenerate combination of the two 
parameters, approximately Qh x with x ~ 1 (i.e. a combination close to the T parameter) 
that is being determined tightly by the elongated ridge of high likelihood. 



flh exp(-ft 6 - /iJ/ 2 ft 6 /ft)(/?,Mpc- 1 ) 
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Figure 3: Contour map of log likelihood in the 0, — h plane for the ACDM (a) and OCDM 
(b) models with n = 1. A[ln(£)] = 1. The most likely values of ft for two given values of 
h, and the corresponding la error bars, are shown. 



We can thus quote stringent constraints on the conditional best value of ft given 
h: ft « (0.7 ± 0.06)/?,- 1 - 2 for ACDM, and ft « (0.75 ± 0.06)h~^° for OCDM (h 50 = 
ffo/SOkms^Mpc- 1 ). 

Table 1 lists the most likely model and it's ln-likelihood (with the zero set arbitrarily) 
for each family of models, within the range ft < 1 and n < 1 and for two values of h. The 
errors are the la conditional uncertainty. We can see from Fig. 2 and the table that, when 
n = 1 is enforced, the best open model is always more likely than the best A model, for 
any given h. 



4.2.2 Tilted Models 

For the tilted ACDM family of models, with or without tensor fluctuations, we use the 
COBE normalization by White and Bunn (1995) (computed only for h = 0.5 and 0.75). 
The results are presented in Figure 4. The best fits and their likelihood are listed in Table 
1. The analysis prefers the highest ft in the range, i.e. , ft = 1, with a nonzero tilt. It also 
prefers the higher value h = 0.75. Nonzero tensor fluctuations increase the likelihood. The 
best fit is obtained at ft = 1, h = 0.75, n = 0.8 and T/S = 7(1 — n). The tight constraint 
is Slh\$n v 0.7 ± 0.08, with v = 1.85 for T/S = 0, and v = 3.8 for T/S = 7(1 - n). 

For fixed ft, this relation can be understood qualitatively as follows: the normalization 
by COBE fixes the amplitude at small wavenumbers, k ~ 0.001, and the velocity data 
constrain the amplitude at k ~ 0.1. The wavenumber corresponding to the peak of the PS 
is proportional to Qh. Therefore, if a good fit is obtained with certain values of h and n, a 
similarly good fit can be obtained with higher h and lower n, or vice versa. The presence 
of tensor fluctuations lowers the amplitude imposed by COBE at small wavenumbers, and 
thus weakens the requirement for a tilt in n. 
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Figure 4: Contours of log likelihood in the Cl — n plane, calculated with h = 0.5 and 0.75 
for TACDM models with and without tensor component. A[ln(£)] = 1. 



The results for the tilted-open family of models (normalized by Sugiyama 1995), for 
the case T/S = only, are presented in Figure 5. The tendency towards large and h is 
similar to the case of tilted-A models, and the maximum likelihood is similar too, but the 
optimal tilt is more pronounced, n = 0.64. The tight constraint in this case is not very 
different either, nh 50 n 177 (0.75 ± 0.08). 

The best-fit power spectra for the families of models discussed above are drawn in 
Figure 6. All the models agree to within ~ 20% for k > 0.1/iMpc -1 , and they differ by 
up to 30 — 50% on larger scales. Figure 7 shows that the most probable model from the 
COBE-normalized CDM models agrees with the T-model (that is COBE independent) on 
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Figure 5: The same as Figure 4 for tilted Open CDM models without tensor component, 
A[\n(C)] = l. 
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Figure 6: The best-fit power spectra for the various CDM models. 

Figure 7: The PS of the most probable COBE-normalized CDM model (solid), and the 
scatter about it following parameter pairs that lie on the innermost likelihood contour 
(Fig. 4, lower-right panel). The COBE-free T model is also marked (heavy solid). The 
PS computed by KD from POTENT density of the same velocity data (independent of 
COBE or models), and their measurement errors, are shown in three bins. 

all scales to within ~ 10%. This is a demonstration of the fact that the velocity data 
contain enough information to constrain the scale of the peak in the power spectrum even 
without the constraint from COBE. Fig. 7 also shows the typical scatter in the PS about 
the best CDM model, using parameter pairs that lie on the innermost likelihood contour 
(Fig. 4, lower-right panel). This scatter is somewhat smaller than in the T model (Figure 
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2b) because of the additional constraint from COBE. 

For comparison, Figure 7 also displays the PS as computed by KD from the POTENT 
smoothed density field recovered from the same Mark III data (independent of COBE or 
models). The results (for all the models tested here) agree within la of the measurement 
errors, and they agree particularly well near k = 0.1 /i Mpc -1 , where the velocity data 
imposes the strongest constraints. The KD spectrum seems to be somewhat steeper. In 
fact it is steeper than any of the CDM spectra discussed here (see Fig. 6), and is roughly 
as steep as the PS predicted by the CHDM model, a 3:7 mixture of hot and cold dark 
matter (KD). The result of the current paper is probably more reliable than KD on large 
scales, because the likelihood method uses all the velocity data including the large-scale 
flows, while the POTENT density field is insensitive to the bulk velocity. On the other 
hand, our results here may be less reliable on small scales because, other than grouping, 
we do not make any correction for nonlinear effects. 

Table 1: CDM Models ft < 1, n < 1 



CDM Model 


ft 


A 


H 


n 


T/S 


ln£ 


Standard 


1 





50 


1 





<o 




1 





40±2 


1 





2.15 


A 


0.70±0.03 


1 - ft 


50 


1 





1.16 




0.415±0.025 


1 - ft 


75 


1 





0.00 


Open 


0.74±0.02 





50 


1 





2.71 




0.54±0.02 





75 


1 





3.66 


Tilted 


1 





50 


0.84±0.04 





2.68 


A or Open 


1 





75 


0.64±0.04 





4.41 


Tilted+tensor 


1 





50 


0.91+0.03 


7(1 - n) 


3.71 




1 





75 


0.80±0.02 


7(1 - n) 


5.31 



5. CONCLUSION AND DISCUSSION 

We have presented a Bayesian method for deriving the power spectrum of mass density 
fluctuations from the Mark III Catalog of Peculiar Velocities. The result is free of galaxy 
"biasing." The method extracts the maximum amount of useful information from the data. 
It is exact to first order, under the assumption of Gaussian fluctuations and Gaussian errors. 
Tests using realistic mock catalogs show that this approximation is adequate because of 
the large-scale coherence of the velocity field and because the large errors that dominate 
on small, nonlinear scales mean that nonlinear effects contribute only weakly to the result. 

Our robust result for the whole family of models examined here is that the PS am- 
plitude at k = 0.1 /i Mpc" 1 is P(k)n 12 = (4.1 ± 0.7) x 10 3 (h^Mpcf, with the PS 
peak in the range 0.03 < k < 0.06 h Mpc -1 , which, for the best-fit model, translates 
to agQ, ' 6 = 0.85 ± 0.1. The errors quoted are crude: they are the typical la uncertainty 
for each of the best-fits within each family of models, and also the typical scatter among 
these best-fit models. Similar results are obtained when using the velocity data alone, 
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as when the additional constraints from COBE are included. Moreover, the results are 
insensitive to the actual choice of model within the general family of CDM models. Note, 
however, that our quoted errors do not include the cosmic scatter that arises from the fact 
that we do not necessarily sample a fair sample of the universe. A crude estimate of the 
cosmic scatter in the PS for the standard CDM model can be found in KD (Fig. 8). It is 
comparable to, and even somewhat larger than, the measurement errors. 

This normalization of the PS is in pleasant agreement with the independent com- 
putation of the PS by KD, which yielded P(k)tt 12 = (4.6 ± 1.4) x 10 3 (h^Mpcf at 
k = 0.1/iMpc -1 , and ergft ' 6 ~ 0.7 — 0.8. Our new results differ at about the 2a level 
from the earlier estimate by Seljak & Bertschinger (1994) based on the earlier Mark II 
sample, of ergft ' 6 = 1.3 ± 0.3. The main improvements since then are that the current 
analysis includes five times denser sampling in a more extended volume, the systematic 
errors such as Malmquist bias are handled better, and a wider span of models is used in 
the likelihood analysis. It may be interesting to note that the current measurement is 
somewhat higher than the completely independent estimate of a similar quantity based on 
cluster abundances, <7 8 ft - 56 ~ 0.57 ± 0.05 (White, Efstathiou, & Frenk 1993), but it is 
only about 2a away. 

The comparisons of the mass ag to the values observed for optical galaxies (~ 0.95) 
and for IRAS 1.2Jy galaxies (~ 0.6 — 0.7) indicate (3 values of order unity to within 25% 
for most galaxy types on these scales (see KD, Fig. 6 and Table 1, for more details). 

A T-shape model, free of COBE normalization, is constrained by the velocity data to 
r = 0.5 ± 0.15. The errors in ag and T are not independent; higher values of T correspond 
to lower normalization, i.e. lower values of ergft ' 6 . This value of V is somewhat higher 
than the canonical values of ~ 0.2 — 0.3 typically obtained from galaxy density surveys 
(e.g. Efstathiou et al. 1992, Peacock & Dodds 1994). 

Within the families of COBE-normalized CDM models, which we have restricted to 
the range ft < 1, the best-fit tilted models with ft = 1 [n = (0.84 ± 0.04)h^ ' ], and the 
best-fit open models with n = 1 [ft = (0.75 ± 0.03)h^ ' ], are found to be more likely than 
the best-fit A models with A = 1 - ft and n = 1 [ft = (0.70 ± 0.03)/?," 1 - 2 ]. 

Our analysis shows that the most likely among all the CDM models has ft = 1, 
h 0.75, and a tilted spectrum of n = 0.8 ± 0.02 with tensor fluctuations of T/S = 
7(1 — n). The most stringent constraints obtained using this likelihood analysis are of 
the sort £lh^ Q n v = 0.7 ± 0.08 with //. = 1.2 for the A models and with v = 3.8 and 1.85 
w/wo tensor fluctuations respectively. For the open models without tensor fluctuations it 
is //. = 1.0 and v = 1.77. 

Our results are consistent with the conclusion of White et al. (1995), who argue for 
tilted CDM based on several data sets including power spectra of galaxy density (Peacock 
& Dodds 1994), cluster correlations, pair-wise velocities and COBE's results. Their best 
fit is ft = 1, h 0.45, n = 0.9, with tensor fluctuations. This model is about la away from 
our best fit (see Fig. 6), but the PS is quite similar; our lower value of n compensates for 
our higher value of h. 

A note of caution about the method. If not enough constraints are imposed, the 
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inversion of Uij (Eq. 5) may be dominated by the noise rather than the signal. For 
example, when we tried to parameterize the PS with a function that did not enforce any 
upper bound on the power on large scales (as is properly enforced in the T model, or when 
COBE constraints are used), then the likelihood analysis preferred unphysically high power 
on large scales. This is, at least in part, a result of noise dominance. An algorithm to 
detect and possibly eliminate this problem is discussed elsewhere (Zaroubi 1996a; 1996b). 
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